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ABSTRACT 

The galaxy group NGC 5044 consists of a luminous giant elliptical galaxy surrounded by a cluster 
of 160 low luminosity and dwarf galaxies, mostly of early type. The cumulative projected radial 
distribution of dwarf galaxies in the NGC 5044 group, unlike distributions of more luminous galaxies 
in rich clusters, does not follow a projected dark matter (NFW) profile. A deficiency or absence of low 
luminosity galaxies is apparent in NGC 5044 within about 350 kpc. Most of the dwarf galaxies in NGC 
5044 entered the virial radius at redshifts z 2 — 3(a//0.25), where a/ = 1/(1 + Zf) is the epoch of 
group formation, and very few entered during redshifts z > 2 — 3(a//0.25). The peculiar, non-NFW 
shape of the projected cumulative dwarf galaxy distribution in NGC 5044 within 350 kpc resembles 
the characteristic cumulative distribution of dark subhalos that are also known to be relatively young. 
Dynamical friction is unlikely to explain the apparent lack of group member galaxies at small radii in 
NGC 5044. 



Subject headings: galaxies: elliptical and lenticular, CD 
feedback 



galaxies: dwarfs - clusters - groups 



NGC 5044 group does not follow the standard projected 
mass distribution of the dark matter halo. Instead, the 
number of dwarf galaxies per unit dark matter appears to 
be lower within a projected radius R ~ 350 kpc. We show 
here that this central shortfall implies that dwarf galaxies 
were underabundant relative to the ambient dark matter 
at redshifts z > 3. This cosmic variation of dwarf galaxy 
density is opposite to that recently proposed to explain 
the increasing number of dwarf galaxies in richer clusters. 

2. OBSERVED GALAXIES IN NGC 5044 

On first inspection, the central elliptical galaxy in the 
NGC 5044 group, with luminosity Lb = 4.5 x 10^° Lb,q, 
completely dominates the optical image of the cluster. A 
closer look reveals a cluster of about 160 low-luminosity 
galaxies surrounding the central elliptical that extends out 
to 500 kpc and beyond (Ferguson & Sandage 1990; FS). 
At 33 Mpc (e.g. Tonry et al. 2001) the absolute magni- 
tudes of the non-central galaxies range from Mb = — 19 
to -13, i.e. 2 X 10^ < Lb < 9 X 10^ Lb,q. About 80 
percent of the non-central group members appear to be 
spheroidal or dE galaxies (FS). Most or all of these galax- 
ies are dwarfs, depending on one's threshold for dwarf des- 
ignation. Unfortunately, cluster membership has not yet 
been fully confirmed with radial velocity measurements, 
but we have removed several Ferguson-Sandage galaxies 
from the membership list - FS68, FS102, and FS137 - 
because their velocities differ by more than ±800 km s~^ 
from that of the central elliptical (NED, Garcia 1993). 
Nevertheless, the remarkably high space density of these 
galaxies and their symmetric distribution largely confirm 
their physical association with the central elliptical (NGC 
5044) . In any case, we assume here that most of the galax- 
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1. INTRODUCTION 

The early evolution of dwarf galaxies is of interest be- 
cause of the possible great antiquity of these galaxies and 
because they lie near the lowest mass scale for hierarchi- 
cal assembly in which low mass galaxies merge to form 
normal galaxies of larger mass. Efforts to understand the 
evolution of these faint galaxies have been frustrated be- 
cause of the difficulty of direct observation at cosmolog- 
ically relevant redshifts. Theoretical models of the for- 
mation of dwarf galaxies are complicated by the poorly 
understood role of supernova feedback and cosmic pho- 
toionization that may limit star formation at certain red- 
shifts. Finally, there is much consternation about N-body 
calculations that predict many more dwarf galaxies than 
are actually observed (Moore et al. 1999). 

We describe here a simple means of determining the 
space density of dwarf galaxies in the early universe by 
observing their projected spatial distribution in galaxy 
groups. Of particular interest are galaxy groups, such as 
NGC 5044, that contain large numbers of dwarf galaxies 
and for which the virial mass is well known from X-ray 
observations. These groups are often dominated by a sin- 
gle, centrally-located luminous elliptical galaxy. Since the 
masses of dwarf galaxies are very small compared to the 
group virial mass, their dynamics are very similar to dark 
matter particles. As dark halos grow by accretion from 
the inside out, so does the spatial distribution of asso- 
ciated dwarf galaxies. Dwarf galaxies currently orbiting 
closest to the central elliptical galaxy must be among the 
oldest galaxies in the group. 

It is therefore of considerable interest that the cumula- 
tive projected radial distribution of dwarf galaxies in the 
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ies regarded as members or likely members by Ferguson & 
Sandage are indeed members. 

Of particular interest is the radial distribution of the 
galaxies in NGC 5044. Carlberg, Yee & Ellingsou (1997) 
and Biviano & Girardi (2003) have shown that the mean 
surface number density of galaxies in richer clusters closely 
follows the NFW profile expected for the underlying dark 
matter (Navarro, Frenk, & White 1997). The luminosities 
of the galaxies considered by Carlberg et al. are much 
larger than those found in NGC 5044. An NFW distri- 
bution for optical galaxies would be expected if galaxies 
entered the virial mass of the cluster evenly over time and 
suffered no significant radial migration by dynamical fric- 
tion not also experienced by the dark matter. 

Remarkably, the projected distribution of the low lu- 
minosity galaxies in NGC 5044 differs dramatically from 
an NFW profile. Instead of the NFW central peak, the 
satellite galaxies in NGC 5044 have a King-like core that 
extends to ~ 150 kpc (Ferguson & Sandage 1990). This 
difference is apparent in Figure 1 where we compare the 
cumulative number of non-central galaxies as a function of 
projected radius R with an NFW mass surface density pro- 
file (Bartelmann 1996). The current virial mass and radius 
of NGC 5044, M^,o w 4 x lO^^ Mq and r^,o = 887 kpc, are 
known from X-ray observations (Buote, Brighenti & Math- 
ews 2004). In Figure 1 both distributions are normalized 
to agree at R = 350 kpc. A Kolmogorov-Smirnov test ap- 
plied within 350 kpc reveals that the observed cumulative 
galaxy distribution in NGC 5044 has a nearly vanishing 
likelihood of being consistent with the NFW profile. The 
normalization that we have chosen, forcing agreement at 
R = 350 kpc, is reasonable because the two distributions 
then agree reasonably well in the range 325 Ss i? ^ 500 
kpc. For R ^ 500 kpc the number of identified galaxies 
falls below the NFW profile. We assume that this dis- 
crepancy is due to an incompleteness in identifying cluster 
members, but the projected radius where the incomplete- 
ness sets in could be somewhat smaller. 

Galaxies that joined the NGC 5044 group at high red- 
shift when the virial radius was smaller have orbits with 
maximum radial excursions near this virial radius. These 
galaxies spend most of their time orbiting near this same 
radius. As we show below, dynamical friction is relatively 
ineffective for dwarf galaxies. Therefore, the numerical 
shortfall of low luminosity galaxies in the core of NGC 
5044 (R < 350 kpc) suggests that proportionally fewer 
galaxies entered the virial radius at high redshift. In the 
following section we develop a simple dynamical model for 
the evolution of orbiting galaxies in NGC 5044 that can 
be used to quantify the relative absence of low-luminosity 
galaxies in the early universe. 

The shallow projected distribution of galaxies around 
NGC 5044 may not be unique. Klypin et al. (1999) 
note the lack of dwarfs in low density regions like the 
Local Group. Smith, Martinez & Graham (2003) find 
that the radial surface density of dwarf galaxies varies as 
^-o.6±o.2 around a sample of isolated elliptical galaxies; 
this is fiatter than NFW except for the innermost regions. 
Madore, Freedman & Bothun (2004) also find that the 
surface density of satellite galaxies near individual ellip- 
ticals declines as while earlier studies by Vader 
& Sandage (1991) and Lorrimer et al (1994) indicate a 



steeper decline, R^^. Dwarf galaxies in the Fornax clus- 
ter have a shallower cumulative radial distribution than 
the giant galaxies (Drinkwater et al. 2001). Very recently, 
Pracy et al. (2004) have discovered a relative deficiency of 
dwarf galaxies in the central regions of the cluster Abell 
2218. The physical scale of the region of dwarf depletion, 
^ 400 kpc, is very similar to that in NGC 5044. In Abell 
2218 the ratio of the projected numbers of dwarf to giant 
galaxies decreases with increasing surface density of giant 
galaxies, as noted by Phillipps et al. (1998). 

3. DYNAMICAL MODEL FOR SATELLITE GALAXIES 

We estimate the expected radial distribution of galaxies 
in NGC 5044 by following the orbits of cluster member 
galaxies in the time-dependent dark matter potential of 
the cluster. The galaxies are regarded as test particles 
without structure. Newly accreted galaxies are assumed 
to arrive at the virial radius having fallen from the local 
turnaround radius. The orbits are not entirely radial since 
galaxies acquire some angular momentum due to attrac- 
tions by neighboring inhomogeneities. After entering the 
virial radius, galaxies orbit until the present time when 
their projected distances from the center of the cluster 
can be compared with the radial distribution observed in 
NGC 5044. 

The first model we consider is a galaxy cluster in which 
the local space density of galaxies is proportional to that 
of the dark matter, as if a fixed fraction of baryons formed 
into low-luminosity galaxies at a very early time. In this 
case the dynamics of dark matter and surviving dwarf 
galaxies are identical. Dark matter halos grow from the in- 
side out by radial accumulation and accretion by mergers. 
The mass assembly histories of dark halos are universal 
and scale invariant, apart from differences in formation 
times. From studies of N-body simulations in a ACDM 
cosmology it has been determined that the mass of viri- 
alized dark matter increases with decreasing redshift z to 
its current value M„_o according to 

= M^,oe-^"f' (1) 

where a/ = l/(l + z/) defines the halo "formation" time 
(Wechsler et al. 2002). We adopt a/ = 0.25 here, but the 
redshifts we derive can be easily scaled to other values of 
aj by keeping afZ constant. At every time or redshift the 
radial profile of dark matter is described by an NFW pro- 
file (Navarro, Frenk & White 1997) that is parameterized 
by the virial mass My and the concentration 

c(M„, z) = 45O(M„/M0)-O-i283(l + z)-' (2) 

(Wechsler et al. 2002), based on a ACDM cosmology with 
no = 0.3, 17a,o = 0.7 and Ho = 70 km s'^ Mpc'^. 

For the idealized circumstance in which dwarf galaxies 
observed in clusters today all formed at some very early 
time, they enter the virialized halo at the same rate as the 
dark matter. With this simple model the virialized mass 
My^i when the ith galaxy first entered the virialized halo 
can be found from 

M„,, 7^M,,o (3) 

where < 7^ < 1 is a random number. The redshift Zi 
at entry can be found from Equation (1) and the corre- 
sponding cosmic time ti{zi) is determined by the ACDM 
cosmology. 
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The virial radius at time U is 



\A-KAiPc,i) 



1/3 



(4) 



where pc,i = SHf/SnG is the critical density and the ap- 
propriate ACDM ovcrdcnsity factor is 

Ai = IStt^ + {n{zi) - 1)[82 - 39(r2(zi) - 1)]. (5) 

The velocity of galaxies arriving at r„_i, 

= (GM„,i/r,,i)l/^ (6) 

is the freefall velocity from the turnaround radius rt — 2r„ , 
and we ignore the small dark energy term, i.e. Ar^ <C 

If the dark matter distribution were perfectly smooth, 
the orbits of dark matter particles and galaxies would 
be perfectly radial. In this case the dark matter halos 
would be more centrally concentrated than the NFW pro- 
file (Bertschinger 1985). However, because of inhomo- 
gcneitics in the dark matter distribution, random velocities 
are induced by mutual gravitational attraction. Conse- 
quently, galaxies arriving at the virial radius have acquired 
some angular momentum that is essential in determining 
the NFW profile. Considering the importance of this ef- 
fect and the large number of available N-body dark matter 
calculations, it is surprising that the angular momentum 
distribution has not been discussed in full detail. We are 
fortunate, however, that Tormen (1997) and Ghigna et al. 
(1998) have plotted the distribution of circularity e = J/Jc 
at r„ where J is the specific angular momentum of dark 
matter particles and Jc = the maximum angular 

momentum at the virial radius; e = 1 for circular orbits 
and e = for radial orbits. 

Based on these results we model the probability density 
for circularity p' (e) = dp/de with a third order polynomial, 
subject to the constraints that p'{0) = p'(l) = and nor- 
malized so the integral of p'{e) from to 1 is unity. The 
resulting distribution is 

p'(e) = ae^-{1.5a+6)e'^ + {0.5a+6)e = Ae^-Be^+Ce (7) 
where 

a = 60(1 - 2(e)) (8) 

and (e) is the mean circularity. The results of Tormen sug- 
gest that (e) « 0.5, for which p'{e) becomes a symmetric 
quadratic. However, the circularity distribution of Ghigna 

et al. is skewed toward larger e so wc also entertain this 
possibility. Random values of e can be found from 
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£ = (A/4)e4 - {B/3)e^ + {C/2)e^. (9) 



All galaxies have velocity Uv,i given by Equation (6) when 
they first enter the virial radius; the angular momentum 
J = JcC is used to allocate the direction of this velocity. 
The galaxy orbits are found by solving 



dr 



= u and 



du 



GM{r) r 



dt at V r 

Within the virial radius the NFW mass distribution is 

J{y) 



(10) 



M{r) = My{t)- 



(11) 



/(c) 

where /(y) = ln(l + y) — ?//(l + y) and y = cr/ry. At 
higher redshifts, when the virial mass grows more rapidly, 
most entering galaxies orbit entirely within the increasing 



virial radius. At lower redshifts, however, some galaxies 
move briefly outside r„(t). In this case we use the exte- 
rior dark matter mass density distribution corresponding 
to the density derived by Barkana (2003) 



P = Pv+ir/ry)' 



-1.3 



r > r„. 



(12) 



The density just beyond the virial radius is found from 

1 dMy 



Pv+ = 



Airrluy dt 



(13) 



Galactic orbits are computed from the time ty^i when 
they first reach the virial radius until the current time 
tn = 13.5 Gyrs. At time t„ the radius of each galaxy n is 
randomly projected on the plane of the sky. 



i^i = [l-(27^-l)2]l/2^ 



(14) 



4. RESULTS 



The thin lines in Figure 1 show the results of two dynam- 
ical models calculated as described in the previous section. 
For these models, and those discussed below, the Monte 
Carlo simulations are continued until 1210 galaxies occupy 
the region R < 350 kpc at time t„. This is ten times the 
observed number in this region. The total number of or- 
bits required to achieve this limit is typically almost twice 
as large. The computed projected radial distributions are 
normalized to agree with the NFW mass distribution at a 
projected radius R — 350 kpc. The light curve (actually 
a histogram) in Figure 1 labeled (e) = 0.5 was computed 
with this constant mean circularity for all orbits during 
the evolution of the dark halo. Since Equation (3) re- 
quires that the arrival of galaxies into the virialized dark 
halo is identical to that of the dark matter, their cumula- 
tive distribution should be identical to the NFW profile. 
The agreement in Figure 1 is seen to be very good indeed, 
particularly considering the simplicity of our model. The 
mean assembly history for the dark halo implicitly allows 
for all past halo-halo mergers. When two halos merge, 
orbiting bound galaxies also merge in a manner that ap- 
proximately conserves the orbital phase space density, and 
this maintains their relative mean locations in the merged 
dark halo potential. 

As the projected radius R approaches the current virial 
radius of NGC 5044, r„,o = 887 kpc, the number of galax- 
ies in the (e) = 0.5 calculation in Figure 1 falls about 10 
percent below the NFW expectation. This deviation prob- 
ably occurs because there are too many (~ 100) galaxies 
at radii r > ry^o at time t„. This discrepancy is insensi- 
tive to the value of the mean circularity (e), i.e. nearly 
identical results are obtained for (e) = 0.3 or 0.7. How- 
ever, if (e) is assumed to increase steadily with the viri- 
alized mass, (e) = My/My^o, the discrepancy is somewhat 
reduced. This result with time- varying (e), illustrated in 
Figure 1 as "var (e)" , has no independent dynamical moti- 
vation. Other small procedural changes, perhaps in Equa- 
tion (6), could be considered to reduce the discrepancy 
near the virial radius, but we have not explored this since 
the results shown in Figure 1 are certainly good enough 
for the region of most interest, R < 500 kpc. 

Figure 2 shows the average halo entry redshift (zj) for 
galaxies in annular bins at time tn for the two dynamical 
models shown in Figure 1 (open and filled circles). The 
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redshift dispersions arc mostly intrinsic, and do not sim- 
ply reflect the finite number of orbiting galaxies in each 
projected bin. As expected, the mean arrival redshift rises 
sharply toward the center of the cluster. Provided the 
space density of low luminosity galaxies is proportional to 
the dark matter density at all redshifts, as imposed by 
Equation (3), the increase in mean redshift with decreas- 
ing projected radius R is robust, essentially independent 
of the value or time variation of (e). In the following we 
assume (e) = 0.5. 

The sense of the discrepancy in Figure 1 between the cu- 
mulative number of galaxies in NGC 5044 and the NFW 
profile indicates that the space density of low luminos- 
ity galaxies relative to dark matter was lower at higher 
redshifts. This discrepancy can be incorporated into our 
models if we relax the assumption of uniform, early galaxy 
formation (Equation 3). We shall consider only models in 
which the space density of dwarf galaxies is constant or 
increases monotonically with time. 

Suppose for example that the satellite galaxies in NGC 
5044 began to form at a redshift Zc when the halo mass 
was Mc and increased as a power law afterward. In this 
case the probability density for arrival of a galaxy into the 
virialized halo is p'{Mv) = Po{My - Mc)«. When p'{My) 
is properly normalized, the virial mass at which a random 
galaxy enters the halo is 

M„,i = Me + (M„,„ - Me)7^Vl+9 (15) 

where 7?. is a random number. This reduces to Equation 
(3) when Mc = 3 = 0. 

Figure 3 shows three evolutionary sequences in which 
galaxies are assumed to begin forming at the earliest pos- 
sible epoch {Mc = 0) and increase in number afterward 
according to a power law, p' oc M^. Of the three possible 
powerlaws, g = 2 is seen to follow best the observed num- 
ber of galaxies in NGC 5044 at small R. As before, all 
solutions are normalized s.t R — 350 kpc. Solutions with 
q = 0.75 (g = 3.5) produce too many (few) galaxies within 
350 kpc, but the q = 2 histogram follows the observed his- 
togram almost perfectly in this region. Nevertheless, none 
of these pure powerlaw models produce the slope change 
at ~ 350 kpc that is apparent in the observations. If the 
q = 2 solution is correct, the observations would have to 
begin to be incomplete just where this slope change oc- 
curs. This is possible but implausible. It seems unlikely 
that Ferguson & Sandage would have missed identifying 
^ 50 member galaxies in the region 350 < i? < 500 kpc. 
The total number of galaxies expected inside the projected 
virial radius 887 kpc at tn is 237, 305 and 352 for orbits 
with q — 0.75, 2 and 3.5 respectively. 

In Figure 4 we plot a better model in which no small 
galaxies form before redshift Zc = 3.32 (Mg = 0.19M^_o), 
followed by a powerlaw rise with q = 0.1. This solution 
fits the observations almost perfectly in i? < 350 kpc and 
continues to increase thereafter in approximate proportion 
to the dark matter. The mean binned redshifts for this dy- 
namical model, plotted as open squares in Figure 2, shows 
that the observations are consistent with very few galaxies 
in NGC 5044 having accretion redshifts > 2. 

Since the powerlaw rise with q = Q.l < 1 is abrupt, 
we also considered a variety of step function models for 
p'{My). In these models g — > and the virial mass at the 



time of galaxy arrival is 

M^,i = Mc + 7^(M„,<, - Mc) (16) 

where the halo mass Mc corresponds to redshift Zc = 
—2ln{Mc/My^o)- The step function solution with Zc = 
2.77 shown in Figure 4 also agrees very well with the NGC 

5044 observations, assuming some small incompleteness at 
R > 350 kpc and few foreground or background interlop- 
ers. 

As a final alternative model we show in Figure 4 a solu- 
tion in which the probability density p'(M„) rises initially 
as a power law until = M* and is constant thereafter. 
For models of this type the mass at which galaxies arrive 
in the halo is given by 

M„,i =M*(7^/7^*)Vp+l (17) 

if 

n<Tl, = [l + {p+ l){My,o - M*)/M,]-^ (18) 

and 

M„,i = M41 - (1 - TZ/TZ,)/{p+l)] (19) 

if 7^ > 7^*. A solution with p = 2.5 and M, = 0.4M„,o 
{z* = 1.83) is shown in Figure 4. This solution also fits 
the NGC 5044 data remarkably well within 350 kpc and 
corresponds to a modest incompleteness in the observed 
number of galaxies at larger radii. 

Overall, we are struck by the ease at which our vari- 
ous dynamical models having a relative absence of high 
redshift dwarf galaxies fit the histogram observed by FS 
within ^ 350 kpc. This excellent agreement gives us con- 
fidence that the FS data are reasonably complete in the 
dense central region of NGC 5044 and that the number of 
dwarf galaxies per unit dark matter is indeed much lower 
for entry redshifts z > 2.5/(a//0.25) than afterwards. 

In addition, we note that the cumulative radial pro- 
jected distributions of dark subhalos in N-body calcula- 
tions (e.g. de Lucia, G. et al. 2004; Gao et al 2004) are 
also shallower than that of the underlying dark halo and 
have histograms very similar to the observed cumulative 
number of dwarfs in NGC 5044. This is reasonable since 
subhalos are continuously destroyed, so surviving subhalos 
are relatively young (Gao et al. 2004). However, the evo- 
lution of dark subhalos and dwarf galaxies differs because 
the subhalos are known to be disrupted and dissolved by 
tidal forces and dynamical friction. 

5. DYNAMICAL FRICTION 

Can the scarcity of high redshift dE galaxies in the NGC 

5044 group be understood as a depletion of central galaxies 
by dynamical friction? The answer is not so obvious since 
dynamical friction increases with the orbiting mass and 
the initial dark mass associated with each dwarf galaxy 
may have been rather large. However, if the central region 
of NGC 5044 is depopulated by dynamical friction, it may 
be repopulated by galaxies that decay by dynamical fric- 
tion from initially larger orbits. To address this question 
we have made a series of orbital calculations including an 
additional dynamical friction term in the second Equation 
(10), 

= -u • 47rlnAG^Mpu-^[ei{{X) - ^Xe"^']. 

(20) 



5 



(Chandrasekhar 1943). Here p{r,t) is the local density 
of halo dark matter, u = |u| is the galactic velocity, 
X = u/2^/^cr and a{r, t) is the mean velocity dispersion of 
dark matter particles in the NGC 5044 halo. According to 
Hoeft, Miickct & Gottlobcr (2004), the dispersion velocity 
of dark matter is approximately isotropic with maximum 
value (Jm = 1-25 ± 0.20(GM„/r„)V2 ^nd radial distribu- 
tion a{r)/am = [2/(^-°-^ + C"-^)]^/^ where ^ = 2y and 
y = cr/r,,. Wc choose InA = 3 as suggested by Zhang et 
al. (2002) and others. Because of the high density and 
compact nature of dwarf galaxies, we approximate the or- 
biting galaxies with point masses. 

In selecting galaxies wc assume random halo entry red- 
shifts given by Equation (3). The mass of the ith galaxy is 
Mi — TbLb^i where is an adjustable mass to hght ra- 
tio and the galactic luminosity is randomly selected from 
the power law part of the Schechter luminosity function 
dN/dLk oc Lg^-^ (which is satisfied by NGC 5044), 

L5, = [L-0-2_7e(L^o-2_i-o.2)^-5^ (21) 

Here Li = 2 x lO"^ Lb,q and L2 = 9.2 x 10^ Lb,q span 
the range of non-central galaxies in NGC 5044. 

The ratio of stellar to dark mass in dwarf galaxies is un- 
certain since mass to light ratios T b can only be observed 
in the central region containing luminous stars. Observa- 
tions of several dE galaxies by Geha et al. (2002) indicate 
3 ^ Tb ^ 6, but the total value could be higher, depend- 
ing on the severity of tidal truncation. Consequently, wc 
computed galactic orbits using the friction term (20) with 
values of Tb = — 100. Circularities were determined by 
Equation (9) assuming (e) = 0.5 . 

For all values of Tb considered the cumulative num- 
ber of galaxies at time t„ closely follows the NEW profile 
(normalized at i? = 350 kpc) and resembles very closely 
the (e) = 0.5 solution in Eigure 1. For Tb = 4 only 
~ 2 percent of the dwarf galaxies merge into the central 
E galaxy, but this increases to 4 or 15 percent if Tb = 10 
or 100 respectively. Even with these very large Tb. the 
loss of galaxies by inward migration and central mergers 
does not appreciably alter their radial distribution at time 
tn- For all Tb considered the redshift distribution of sur- 
viving galaxies at i„ is essentially the same as that shown 
in Figure 2 for the (e) = 0.5 solution without dynamical 
friction. 

However, when Tb <; 30 we find that the average lu- 
minosity {Lb) of surviving galaxies decreases. Since the 
drag force is proportional to Mdu/dt oc oc (Tb^b)^, 
the onset of dynamical friction is rather abrupt as the 
mass OT Lb of orbiting galaxies increases. This may ex- 
plain why the luminosity distribution of low-luminosity 
galaxies observed in NGC 5044 cuts off rather abruptly at 
Lb ~ L2. In orbital calculations using Equation (3) and 
including dynamical friction, the average luminosity of sur- 
viving satellite galaxies (Lb) drops by 10 or 20 percent if 
Tb = 15 or 30 respectively. For Tb = 65, (Lb) is only half 
as large as the average galaxy observed in Li < Lb < L2 
because a larger fraction of these more massive galaxies 
have computed orbits that merge at the center under dy- 
namical friction. We conclude that the mass to light ratios 
of the (possibly tidally truncated) low luminosity galaxies 
in NGC 5044 cannot be much larger than ~ 30. Summa- 
rizing, (1) for reasonable dwarf galaxy mass to light ratios 
Tb dynamical friction is not very important, and (2) the 



main effect of increasing T b above ~ 30 is to deplete the 
more luminous dwarfs relative to those observed in NGC 
5044. 

6. DISCUSSION AND CONCLUSIONS 

We emphasize the distinction between the evolution of 
group member galaxies and dark matter subhalos, which 
has been a source of confusion in recent discussions. Sub- 
halos arc formed by dark matter inhomogcneitics that arc 
accreted into the main halo. As these subhalos undergo 
tidal disruption and dynamical friction, they melt into 
the main cluster halo so surviving subhalos at any time 
are comparatively young. However, the dense baryonic 
galactic cores in siibhalos are very deeply bound and arc 
much more resistant to dynamical disruption than the sur- 
rounding subhalos. The dynamical distinction between 
ephemeral subhalos and long-lived galactic cores has been 
discussed for massive clusters by Gao et al. (2004) who use 
N-body and semi-analytic methods to reproduce both the 
shallow, non-NEW cumulative distribution of dark subha- 
los and the nearly NFW distribution of surviving galaxies 
(Carlberg, Yee & EUingson 1997). Within the resolution 
limit of the calculation of Gao et al., the surviving galax- 
ies typically have lost their extended dark halos by tidal 
disruption and this is likely to be true in NGC 5044 as 
well. 

We also emphasize that the final galactic distributions 
calculated here refer only to cluster member galaxies that 
survive today. Many more massive galaxies that did not 
survive (and are therefore not included in most orbits con- 
sidered here) merged to form the giant elliptical at the 
center of NGC 5044. An integral over a Schechter lumi- 
nosity distribution from the most luminous orbiting galaxy 
in NGC 5044 Lb = L2 to infinite Lb is comparable to the 
luminosity of the central galaxy (Jones, Ponman & Forbes 
2000). This is consistent with the notion that the cen- 
tral elliptical was formed by mergers of all first generation 
galaxies with Lb > L2. 

The relative absence of dwarf galaxies that entered the 
virial radius of NGC 5044 at high redshifts is exactly oppo- 
site to the behavior that has been proposed to explain the 
excess of dwarf galaxies in rich clusters (Klypin et al. 1999; 
Bullock, Kravsov & Weinberg 2000; Trentham, Tully & 
Verheijen 2001; Tully et al. 2002; Tassis et al. 2003). Tully 
et al. (2002) argue, for example, that the relative lack of 
dwarf galaxies in the Ursa Major cluster (M„ = 4 x 10^^ 
Mq) relative to the Virgo cluster (M„ = 8 x lO^^ Mq) 
can be understood if most of the dwarfs in Ursa Major 
formed after cosmic reionization at redshift Zion ^ 6 when 
the heated baryonic gas was less able to form stars. Late 
forming dwarfs also have lower dark matter (and gas) den- 
sities, further suppressing star formation. Conversely, it is 
argued that the Virgo cluster is much older and its (more 
numerous) dwarfs formed before reionization. 

Evidently, dwarf galaxy formation and the faint end of 
the galaxy luminosity function may vary with space and 
time due to environmental influences. The relatively small 
number of surviving dwarf galaxies that entered the virial 
radius of NGC 5044 during redshifts 6 > z > 2 — 3 could be 
cixplained if stars were unable to form in small galaxies un- 
til z ~ 2— 3. If the dark matter potential in low mass halos 
is sufficiently shallow, heating by early supernova feedback 
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or cosmic ionizing radiation may have kept the baryonic 
gas from coohng into stars. If no gas accumulated in or 
near low mass dark halos prior to redshift .2 ~ 2 — 3, only 
the dark halos of dwarf galaxies entered the virial radius 
of NGC 5044. To explain the existence of dwarf galaxies 
in NGC 5044 at smaller redshifts, it must be assumed that 
some baryonic gas eventually accreted into these halos and 
formed into the stars observed today. This scenario is pos- 
sible, but unlikely. According to Susa & Umemura (2004) , 
star formation in dwarfs with viral masses greater than 
~ 10^ Mq is unaffected by cosmic UV heating. For typi- 
cal stellar mass to light ratios Tb,* ~ 5, the stellar masses 
of the NGC 5044 dwarfs all exceed ~ 10^ Mq and their 
original virial masses were ~ 10 times larger. Therefore 
star formation in the dwarfs of NGC 5044 may be impaired 
but not completely shut down by cosmic UV heating. This 
is consistent with observations of nearby dwarf galaxies, 
many of them less massive than those in NGC 5044, that 
show no evidence that the star formation rate changed ei- 
ther at the epoch of reionization or at redshifts of ~ 2 — 3 
(Grebel & Gallagher 2004) 

The apparent deficiency of high redshift dwarf galax- 
ies in NGC 5044 can also be understood if these galaxies 
retained most of their gas at all relevant redshifts, pro- 
vided very few baryons in these galaxies formed into stars 
before z ~ 2 — 3, as implied by the slow and sustained 
star formation in many well-observed local dwarf galaxies 
(e.g. Grebel & Gallagher 2004). Star formation in gas- rich 
dwarfs could be reduced either by proximity to NCC 5044 
or by gas removal after the dwarf entered the virial radius 
of NGC 5044. Suppose for example that some disturbing 
process occurred during formation of the central galaxy 
(NGC 5044) such as UV ouiissioii, outflowing winds, 
AGN activity etc. - that suppressed star formation in 
nearby dwarf galaxies. A necessary condition for this sce- 
nario is that the central formation starburst of NGC 5044 
preceded the formation of most of the stars in the sur- 
rounding system of dwarf galaxies. As discussed already, 
even the least massive of the dwarfs identified in NGC 5044 
is massive enough to retain ionized gas at ~ 10* K, but 
gas heated to this temperature in gas-rich dwarf galaxies 
may expand sufficiently in the dwarf potential to greatly 
increase the efficiency of ram-pressure stripping after the 
galaxies entered the virialized hot gas in the NGC 5044 
group. In this case the stellar component of high redshift 
dwarfs in NGC 5044 would be optially faint. One possible 
difficulty with this scenario is that the redshift range in 
Figure 2 corresponding to the reduced numbers of dwarfs, 
2 — 3^2^6, corresponds to almost 2 Gyrs and it is 
unlikely that an ionizing starburst (or AGN activity) can 
be maintained for this length of time. If the gas in dwarf 
galaxies is cooler than ~ 10** K, it is unclear if it can be 
ram-stripped by the hot gas in NGC 5044. 

Alternatively, could two-body interactions with larger 
galaxies, dark galaxies or subhalos at early times have 
scattered the orbits of dwarf galaxies originally near the 
center to larger radii in the NGC 5044 group? Probably 
not. Normal two-body galactic interactions are strongly 
masked by the large background dark matter density in 
NGC 5044 that dominates the force on each dwarf galaxy. 
For example, the concentrated baryonic mass of the cen- 
tral massive galaxy in this group is overwhelmed by the 



ambient NFW dark mass beyond a radius of only 20 kpc. 
If galaxies are inefficient orbital scatterers, perhaps the 
softer potentials of massive dark subhalos within r ~ 350 
kpc in NGC 5044 strongly perturbed the orbits of dwarfs 
near the center. Unfortunately, subhalo perturbations are 
not likely to be effective since (1) X-ray observations and 
detailed N-body calculations verify that the dark matter 
orbiting in the central region of NGC 5044 is not selec- 
tivly depleted in this way and (2) normal galaxies in the 
corresponding central regions of rich clusters are not scat- 
tered outward by dark subhalos. This is verified by the 
galaxy plus subhalo simulations of Gao et al. (2004) and 
by direct observation (Carlberg et al. 1997). The NFW 
distribution of luminous galaxies in rich clusters can be 
understood if these galaxies formed well in advance of the 
cluster potential. 

We appeal to observers to study NGC 5044 and other 
similar groups in more detail. Group membership, partic- 
ularly at large radii, must be confirmed with accurate ra- 
dial velocities. Measured orbital velocities must also sup- 
port the relative absence of dwarf galaxies at small radii. 
Group memberships should be determined out to the virial 
radius and beyond. If our interpretation of the incom- 
pleteness of the Ferguson-Sandage data is correct, another 
50-100 members in NGC 5044 remain to be discovered in 
400 ^ R <> 900 kpc (Figure 1). If possible, estimates of 
the stellar ages in dwarfs as a function of projected radius 
could constrain our argument for a relative deficiency of 
old dwarfs having accretion redshifts 6 ^ z ^ 2 — 3. An 
extension of group membership to dwarf galaxies of even 
lower luminosity could also be useful. 

Our main conclusions are: (1) The cumulative pro- 
jected radial distribution of dwarf galaxies in the NGC 
5044 group does not follow a projected NFW profile, un- 
like observed distributions of more luminous galaxies in 
rich clusters. (2) Inside about 350 kpc in NGC 5044 there 
is a deficiency or absence of low luminosity galaxies rel- 
ative to the NFW profile. (3) Most of the dwarf galax- 
ies in this inner region of NGC 5044 that would have en- 
tered the virial radius at redshifts 2;/(a//0.25) > 2 — 3 
are largely absent; most of the surviving dwarf galax- 
ies in NGC 5044 first passed through the virial radius 
at redshifts z < 2 — 3(a//0.25), where a/ = 1/(1 -|- Zf) 
is the epoch of group formation. (4) The peculiar, non- 
NFW shape of the depleted projected cumulative dwarf 
galaxy distribution in NGC 5044 within 350 kpc resembles 
the characteristic cumulative distribution of dark subhalos 
that are also known to be relatively young. (5) Dynamical 
friction is an unlikely explanation for the apparent lack of 
group member galaxies at small radii in NGC 5044. (6) 
The mass to light ratio of dwarf galaxies in NGC 5044 
is unlikely to exceed ~ 30 since otherwise the more lu- 
minous observed dwarf galaxies would have been lost by 
dynamical friction. (7) Our finding that dE galaxies did 
not appear until moderate redshifts is qualitatively consis- 
tent with the low mean stellar age of ^ 5 Gyrs observed in 
dE galaxies (Geha et al 2002), i. e. many or most dwarf 
galaxies are younger than typical giant elliptical galaxies. 



Studies of the evolution of hot gas in elliptical galaxies 
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Fig. 1. — The irregular heavy line curve labeled "NGC 5044" is a cumulative histogram of the projected radii of low luminosity galaxies 
in the galaxy group NGC 5044. The heavy dashed line is the cumulative projected mass of an NFW profile normalized to agree with the 
NGC 5044 histogram at R = 350 kpc. The two thin lines are dynamical models (actually histograms) of the cumulative radial distribution of 
dwarf galaxies in NGC 5044 if the galaxies entered the virialized hales in proportion to the dark matter. The circularity is constant for the 
(e) = 0.5 curve and increases with the virial mass from to 1 in the curve labeled "var (e)". 
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Fig. 2. — The mean rcdshift and rcdshift dispersion calculated for member galaxies in NGC 5044 in ammlar bins of projected radius. The 
open circles correspond to the (e) = 0.5 model in Figure 4. The filled circles and open squares correspond respectively to redshifts of model 
galaxies in models "var (e)" (Figure 4) and "zc = 3.32, q = 0.4" (Figure 4). All data points refer to the bins for the open circles but the filled 
circles and open squares have been shifted to the right for clarity. 
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Fig. 3. — The observed histogram for NGC 5044 galaxies and the normalized NFW profile are identical to those in Figure 1. The three 
light lines are histograms of the cumulative number of galaxies in NGC 5044 assuming that they entered the virialized halo with probability 
density p'(M^) oc (M^)« with q = 0.75, 2 and 3.5. 
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Fig. 4. — Three additional models for the entry and dynamics of dwarf galaxies in the NGC 5044 group. For the histogram labeled 
"2:c = 3.32 q = 0.1" no galaxies entered the halo of the group before redshift 3.32 when the virial mass was Mc = O.lQM^.^o, but afterwards 
galaxies entered in proportion to {M^, — Mc)" "'^. The "zc = 2.77 step" curve shows the distribution of galaxies resulting from a step function 
at this redshift. In the galaxy distribution labeled "z, = 1.83 p = 2.5" galaxies entered the halo of NGC 5044 in as M^ '' until redshift 
z* = 1.83 and in proportion to the dark matter afterward. The observed histogram for NGC 5044 galaxies and the normalized NFW profile 
are identical to those in Figure 1. 



